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The heavy quarks constitutes a unique probe of the quark-gluon plasma properties. Both at 
RHIC and LHC energies a puzzling relation between the nuclear modification factor Raa(pt ) and 
the elliptic flow V2(pr) has been observed which challenged all the existing models, especially for 
D mesons. We discuss how the temperature dependence of the heavy quark drag coefficient is 
responsible for a large part of such a puzzle. In particular, we have considered four different models 
to evaluate the temperature dependence of drag and diffusion coefficients propagating through a 
quark gluon plasma (QGP). All the four different models are set to reproduce the same Raa{pt) 
observed in experiments at RHIC and LHC energy. We point out that for the same Raa{pt) one 
can generate 2-3 times more V 2 depending on the temperature dependence of the heavy quark drag 
coefficient. A non-decreasing drag coefficient as T —¥ T c is a major ingredient for a simultaneous 
description of Raa{pt) and V2{pr)- 
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The ongoing nuclear collision programs at Rel¬ 
ativistic Heavy Ion Collider (RHIC) and Large 
Hadron Collider (LHC) energies have created a 
medium that behaves like a nearly perfect fluid. The 
bulk properties of such a matter, called Quark Gluon 
Plasma (QGP), are governed by the light quarks and 
gluons [J, 0- T° characterize the QGP, penetrat¬ 
ing and well calibrated probes are essential. In this 
context, the heavy quarks (HQs), mainly charm and 
bottom quarks, play a vital role since they do not 
constitute the bulk part of the matter owing to their 
larger mass compared to the temperature created in 
ultra-relativistic heavy-ion collisions (uRHIC’s) [3|. 

There are presently two main observables related 
to heavy quarks that have been measured at both 
RHIC and LHC energies. The first one is the so- 
called nuclear suppression factor Raa that is the ra¬ 
tio between the px spectra of heavy flavored hadrons 
(D and B) produced in nucleus + nucleus collisions 
with respect to those produced in proton + pro¬ 
ton collisions. More specifically at RHIC until re¬ 
cently has not been possible to measure directly D 
and B but only the leptons through their semilep- 
tonic decays. The other key observable is the elliptic 
flow V 2 = ( cos(2ip p )), a measure of the anisotropy 
in the angular distribution that corresponds to the 
anisotropic emission of particles with respect to the 
azimuthal angle 4> p . Despite their large mass, exper¬ 
imentally measured nuclear suppression factor Raa 
and elliptic flow V 2 of the heavy mesons are com¬ 
parable to that of light hadrons EMI. This is in 
contrast to the expectations drawn initially from the 
perturbative interaction of HQs with the medium 
which predicted a Raa ~ 0.6 for charm quarks, 
Raa ~ 0.8 — 0.9 for bottom quarks in the central 
collisions EMI at intermediate px- Also the V 2 


was predicted to be much smaller with respect to 
the light hadron ones [ 10 . 

Several theoretical efforts have been made in or¬ 
der to calculate the experimentally observed Raa 
EMI and v 2 E3 f° r the non-photonic single elec¬ 
tron spectra within the Fokker-Planck approach 0- 
HoL [Ta . [ 2 TI fril IHHl [3514381 and relativistic Boltz¬ 
mann transport approach [n| [291433 . [H, [HU . Fur¬ 
thermore, also in a pQCD framework supplemented 
by the hard thermal loop scheme several advances 
have been made to evaluate realistic Debye mass 
and running coupling constants 16, 26 and also 
three-body scattering effects ES I2ll vPA IHll have 
been implemented . It has been show in |39( that 
the inclusion of both elastic and inelastic collisions 
within a dynamical energy loss formalism reduces 
the gap between the theoretical and exp erimental 
results for Raa as px > 5 — 10 GeV [ifl |4l|. Sev¬ 
eral other improvements have been prop osed to ad¬ 
vance the description of the data [43-[44J] . Inter¬ 
action from AdS/CFT [60j have also been imple¬ 
mented [24], El [HU to study the heavy flavor dy¬ 
namics at RHIC and LHC. Essentially all the mod¬ 
els show some difficulties to describe simultaneously 
both Raa{pt) and V 2 (jpx) and such a trait is not 
only present at RHIC energy but also in the results 
coming from collisions at LHC energy (20j |. 

In this letter we will address the impact of the 
temperature dependence of the interaction (drag co¬ 
efficient) on both Raa and V 2 relation. For this we 
are considering four different models having different 
T dependent drag coefficients. For the momentum 
evolution of the HQ, we are using 3+1 D Langevin 
dynamics. We notice that the several approaches 
and modelings of the HQ in-medium interaction dif¬ 
fers significantly for the T dependence of the drag 
coefficient they entail. One can go from a T 2 depen- 
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dence of the AdS/CFT approach to a drag coeffi¬ 
cient that even increase with decrease T. The aim of 
this letter is to show that, while generally a smaller 
Raa(pt ) corresponds to larger V 2 (pr) ; the specific 
T dependence of the drag can strongly modify such 
an amount of v?( pt), even if the models are tuned to 
reproduced the same Raa(pt) observed experimen¬ 
tally. Our analysis shows that it is quite unlike that 
a drag with T 2 dependence can generate larger ellip¬ 
tic flow as the one observed experimentally at both 
RHIC and LHC. Instead a nearly constant drag or 
an increasing one as T —> T c strongly quenches the 
puzzling Raa{pt) — V 2 (pt) relation. 

The standard approach to HQ dynamics in the 
QGP is to follow their evolution by means of a 
Fokker-Plank equation solved stochastically by the 
Langevin equations. The relativistic Langevin equa¬ 
tions of motion for the evolution of the momentum 
and position of the heavy quarks can be written in 
the form 

dxi = ^ dt , 

jb 

dpi = — Tpidt + CijPjVdt (1) 

where dxi and dpi are the shift of the coordinate 
and momentum in each time step dt. T and Cij are 
the drag force and the covariance matrix in terms 
of independent Gaussian-normal distributed random 
variables p,P(p) = (27r) _3//2 e _p / 2 , which obey the 
relations < piPj >= 8ij and < pi >= 0, respectively. 
The covariance matrix is related to the diffusion ten¬ 
sor, 

Cij = y/2B Q P± + VmPl, ( 2 ) 

where PA = Sij -PiPj/p 2 and P} g = PiPj/p 2 are 
the transverse and longitudinal projector operators 
respectively. Under the assump tion, Bq = B\ = D , 
Eq (0 becomes Cij = i/2 D(p)Sij . Such an assump¬ 
tion strictly valid only for p —> 0 , is usually employed 
at finite p in a ppl ication for heavy quark dynamics 

in the QGP EE US HI . 

We will discuss our results in terms of the drag co¬ 
efficient r, but we remind that it is related to the dif¬ 
fusion coefficient by the fluctuation-dissipation the¬ 
orem that within a Langevin approach reads D = 
VET , for the case of the post-point Ito realization 
of the stochastic integral[45j. In the post-point dis¬ 
cretization the diffusion coefficients have to be used 
at the momentum argument p+dp, where dp is the 
increment from a pre-point Ito (Euler) time-step ac¬ 
cording to Eq. 0 . 

The solution of the stochastic Langevin equation 
needs a background medium describing the evolu¬ 
tion of the bulk QGP matter. To describe the the 
expansion and cooling of the bulk and its elliptic 
flow V 2 (pt) at both RHIC and LHC, we have em¬ 
ployed a relativistic transport code with an initial 


condition given by a standard Glauber model and 
with an evolution at fixed v/s = 0.16 (similarly to 
viscous hydro) see Ref.s FfOl M0l | for more details. 

Our objective is to demonstrate the effect of the 
temperature dependent interaction (drag coefficient) 
on the Raa and V2 obtained from different models. 
More specifically we investigate at fixed Raa how 
the V 2 is built up under various temperature depen¬ 
dence of the interaction. For this purpose we con¬ 
sider four different modelings to calculate the drag 
and diffusion coefficients which are the key ingredi¬ 
ents to solve the Langevin equation. Such models 
have to be considered merely as an expedient-device 
to generate different T dependence of the T(T) but 
the results and conclusions deduced will be much 
more general because they do not depend on the way 
the r(T) has been obtained. In this sense within a 
Fokker-Planck approach it is not relevant if the drag 
and diffusion coefficients has been evaluated consid¬ 
ering only collisional or radiative loss. 

Model-I (pQCD) - The elastic interaction of heavy 
quarks with the light quarks, anti-quarks and gluons 
in the bulk has been considered within the frame¬ 
work of pQCD to calculate the drag and diffusion 
coefficients. The scattering matrix A4 gH Q y M. q HQ 
and M.qHQ are the well known Combridge matrix 
that includes s,t,u channel and their interferences 
terms 0. The divergence associated with the t- 
channel diagrams due to massless intermediate par¬ 
ticle exchange has been shielded introducing the De¬ 
bye screening mass mn = \fAnoT s T. The tempera¬ 
ture dependence of the coupling [581 ]: 


g~ 2 (T) = 2/ 0 ln 


2t tT 
aT c 


^In 


In 


2ttT 

aT r 


(3) 


where /3o = (11 — 2TV//3)/167 t 2 , fd\ = (102 — 
38iV//3)/(167r 2 ) 2 and a = 1.3. Nf is the number 
of flavor and Tc is the transition temperature. 

Model-II (AdS/CFT) - We have also considered 
the drag force from the gauge/string duality [59j |. 
namely the conjectured equivalence between con¬ 
formal N=4 SYM gauge theory and gravitational 
theory in Anti de Sitter space-time i.e. AdS/CFT. 
By matching the energy density of QCD and SYM, 
which leads to Tsym = Tqcd/ 3 j , and the string 
prediction for quark-antiquark potential with lattice 
gauge theory which gives 3.5 < A < 8 [60]. One 
finds; 


conf 


= c- 


T; 


QCD 

M r 


(4) 


where C = = 2.1 ± 0.5. The correspond¬ 

ing diffusion constant D can be obtained from the 
fluctuation-dissipation relation. Studies of heavy fla¬ 
vor momentum evolution within the Langevin dy¬ 
namics using AdS /CFT can be found in Ref. [IE Eg ■ 
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Model-Ill (QPM) - The third model recently ap¬ 
plied to estimate the heavy flavor transport co¬ 
efficients is inspired by the quasi-particle model 
(QPM) IM3. The QPM approach is a way to 
account for the non-perturbative dynamics by T- 
dependent quasi-particle masses, m q = 1/3 g 2 T 2 , 
m g = 3/4< 7 2 T 2 , plus a T-dependence background 
field known as bag constant. Such an approach is 
able to successfully reproduce the thermodynamics 
of 1QCD [63] by fitting the coupling g(T). To eval¬ 
uate the drag and diffusion coefficients we have em¬ 
ployed QPM tuned to the thermodynamics of the 
lattice QCD [64j . Such a fit lead to the following 
coupling [631 ]: 


9 2 {T) 


48tt 2 

[(117V C - 2N f )ln[\(% - ^)] 2 


(5) 


where A=2.6 and T/T s = 0.57. 

Model-IV (otQPM(T),m q = m g = 0/ - To have a 
different set of drag and diffusion coefficients we are 
considering a case where the light quarks and glu¬ 
ons are massless but the coupling is from the QPM 
which obtained from the fit to the lattice data. This 
case has to be mainly considered as an expedient to 
have a drag which decreasing with T as obtained for 
example in the T-matrix approach [3i. fl5i. IflisS . 



FIG. J| Variation of drag coefficient with respect to tem¬ 
perature. 


In the following, except for the case of AdS/CFT, 
we have calculated the drag coefficient numerically 
from the scattering matrix of the model by means of 
the standard definition of drag @, see also Ref. [67] 
for a recent detailed description of the calculation of 
the transport coefficients for heavy quarks. 

The variation of the drag coefficient with re¬ 
spect to temperature at p=100 MeV obtained within 
the four different models discussed above has been 
shown in Fig[T] The behaviors remain quite similar 


also at high momentum but with different magni¬ 
tude. These rescaled drag coefficients can repro¬ 
duce almost the same Raa at RHIC energy. In 
AdS/CFT case the drag coefficient is proportional 
to T 2 whereas in otQpM{T),m q = m g = 0 case the 
drag coefficient decrease with T due to the strong 
coupling at low temperature. It may be mentioned 
here that the drag coefficient obtained from the T- 
matrix stum is almost constant or slightly de¬ 
creasing with temperature. 

We mention that the drag coefficient increases 
with temperature when the system behaves like a 
gas. For a molecular liquid the drag coefficient de¬ 
creases with increasing temperature (except in a 
very few cases) because a significant part of the ther¬ 
mal energy goes into making the attraction between 
the interacting particles weaker, allowing them to 
move more freely and hence reducing the drag coef¬ 
ficient. The drag force of the partonic medium with 
non-perturbative effects may decrease with increas¬ 
ing temperature as shown in Ref. UH, @ because 
in this case the medium interacts strongly more like 
a liquid. 

In order to study the impact of the temperature 
dependence of the drag coefficient presented in the 
previous sections on the experimental observables, 
we have calculated the nuclear suppression factor, 
Raa, using our initial charm and bottom quark dis¬ 
tributions at initial time t = t, and final time t = tj 

at the freeze-out temperature as Raa{p) = J§rrj- 
Along with Raa we evaluate the anisotropic momen¬ 
tum distribution induced by the spatial anisotropy 
of the bulk medium and defined as 



( 6 ) 


which measures the momentum space anisotropy. 

We have performed simulation of Au + Au col¬ 
lisions at Va = 200 AGeV for the minimum bias 
using a 3+1D transport approach HE S3, H3]. The 
initial conditions for the bulk evolution in the coordi¬ 
nate space are given by the Glauber model condition, 
while in the momentum space we use a Boltzmann- 
Juttner distribution function up to a transverse mo¬ 
mentum pt = 2 GeV and at larger momenta mini-jet 
distributions as calculated within pQCD at NLO or¬ 
der [23|. At RHIC energy, Au + Au at y/s = 200, 
the maximum initial temperature of the fireball in 
the center is T) = 340 MeV and the initial time 
for the fireball simulations is t* = 0.6 fm/c (accord¬ 
ing to the criteria t» • Tj ~ 1). The heavy quarks 
distribution in momentum space are distributed in 
accordance with the charm distribution in pp colli¬ 
sions that have been taken from ref Si where in 
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FIG. 2: Comparison of the nuclear suppression factor, 
Raa, as a function of pr, obtained within the Langevin 
(LV) evolution for the four different cases, with the ex¬ 
perimental data at RHIC energy. 


FIG. 4: Comparison of the nuclear suppression factor 
Raa vs V2, obtained within the Langevin (LV) evolution 
for the four different cases, with the experimental data 
at RHIC energy at pt = 1.3 GeV. 



FIG. 3: Comparison of the elliptic flow, V 2 , as a function 
of pr, obtained within the Langevin (LV) evolution for 
the four different cases, with the experimental data at 
RHIC energy. 


the coordinate space they are distributed according 
to N coU . 

The solution of the Langevin equation has been 
convoluted with the fragmentation functions of the 
heavy quarks at the quark-hadron transition tem¬ 
perature T c to obtain the momentum distribution of 
the D and B mesons. For the fragmentation, we use 
the Peterson fragmentation function: 


where e c = 0.04 for charm quarks and e c = 0.005 for 
bottom quark. 


In Fig. [2] we have plotted Raa as a function of 
Pt for the four different cases obtained within the 
Langevin dynamics at RHIC energy. As we men¬ 
tioned, we try to reproduce the same Raa in all the 
cases by rescaling the drag and diffusion coefficients. 
We remind that RHIC data and calculations refer to 
the single electrons from the semileptonic decay of 
D and B mesons. The V 2 for the same Raa has 
been displayed in Fig [3] for all cases as a function of 
Pt- Our main striking point is that even if the Raa 
is very similar for all the four different cases, the V 2 
built up is quite different depending on the tempera¬ 
ture dependence of the drag coefficients (see Fig[T]). 
This is because the Raa is more sensitive to the 
early stage of the evolution whereas the V 2 is more 
sensitive to the later stage of the evolution (near T c ). 
Some studies in this direction have been done also 
in the light flavor sector as shown in Ref. |sHz3 
and very recently related to the presence of mag¬ 
netic monopoles [72j]. The larger drag coefficient is 
at low temperature the larger is the V 2 even for the 
same Raa ■ For example in the region of the peak for 
V 2 (pr) we see a difference of about a factor 2.5 going 
from a T 2 dependence, like AdS/CFT to a inverse 
T dependence as it can occur in a liquid. This last 
case or at least a nearly constant drag appears to be 
very much favored by the comparison with the data. 

This study suggests that the correct temperature 
dependence of drag coefficient has a crucial role for 
a simultaneous reproduction of Raa and V 2 - The 
reason for such a relation between the two observ¬ 
ables is that a small Raa (strong suppression) can be 
generate very quickly at the beginning of the QGP 
lifetime, i.e. at high T. However such a strong inter¬ 
action will not be accompanied by a build-up of V 2 
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because the bulk medium has not yet developed a 
sizeable part of its elliptic flow. On the contrary to 
generate a large v-i one needs that there is a strong 
interaction with the medium at later stages of the 
QGP lifetime in order to match the build-up of both 
Raa and V 2 - The experimental data seem to clearly 
suggest that the drag of the medium cannot decrease 
with large power of T otherwise the interaction will 
be relatively weak just when a strong interaction 
would make possible the build of the anisotropy in 
momentum space. It can be here mentioned that 
the drag coefficient is almost constant with res pec t 
to temperature in the T-matrix case f EE HE HU • 
However also a QPM can be considered quite close to 
the data given that we have not included the coales¬ 
cence mechanism that would shift the V 2 (pr) in all 
the cases considered by about a 20-25% upward. In 
Fig 0] we have introduced a new plot Raa vs V 2 at a 
given momentum (pr= 1-3 GeV) to promote the im¬ 
portance of simultaneous reproduction of Raa and 
V 2 - Fig [4] highlights how the V 2 (pr) built up can 
differ up to a factor of around 2.5 (in the region 
of peak), for the same Raa(pt ), depending on the 
temperature dependence of the drag coefficient. 



FIG. 5: Comparison of the nuclear suppression factor, 
Raa , as a function of pr, obtained within the Langevin 
(LV) evolution for the four different cases, with the ex¬ 
perimental data at LHC energy. 


We have also extended our calculation to study 
Raa and V 2 at LHC performing simulations of 
Pb + Pb at y/s = 2.76 ATeV energy. In this case 
the initial maximum temperature in the center of 
the fireball is T 0 = 510 MeV and the initial time for 
the simulations is r 0 ~ 1/T 0 = 0.3 fm/c. In Fig.[5]we 
show the Raa as a function of p t for the four differ¬ 
ent cases obtained within the Langevin dynamics at 
LHC energy. As we mentioned, we reproduce similar 
Raa in all the cases by rescaling the drag and dif¬ 
fusion coefficients. The elliptic flow V 2 for the same 



FIG. 6: Comparison of the elliptic flow, V 2 , as a function 
of pr, obtained within the Langevin (LV) evolution for 
the four different cases, with the experimental data at 
LHC energy. 


Raa has been plotted in Fig [6] for all cases as a 
function of pt ■ Similarly to the RHIC case we get 
a similar trend for the Raa vs V 2 depending on the 
T dependence drag coefficients. 

However, as pointed out in Ref (34j, for charm 
quarks, which have a moderate M/T ratio, a sig¬ 
nificant deviation with respect to the Brownian 
Langevin dynamics can be expected. In this case the 
full solution of the Boltzmann integral i.e. without 
the assumption of small collisional exchanged mo¬ 
menta, leads in general to a large V 2 (pr)- Such an 
effect depends on the anisotropy of the microscopic 
scattering and can not be studied in term of only the 
drag coefficient. It is however an effect that in gen¬ 
eral can be expected to be of the order of about 20% 
and does not modify the systematic studied here. A 
further effect that is involved in the study of HQ ob¬ 
servable is related to the hadronization process. If 
the possibility of the coalescence process is included 
there is a further enhancement of the V 2 (pr) of about 
a 20 — 25% H, [Til [73j. Also the hadronic rescatter¬ 
ing may play a role in enhancing the V 2 (pr) without 
modifying the Raa{pt) fH). This however would 
generate a similar shift for all the cases discussed 
hence not affecting the discussed pattern entailed 
by r(T). The impact of Boltzmann dynamics and 
hadronization by coalescence are larger at LHC and 
can be lead to a better agreements with the data for 
the case Q;qpm(T) and QPM but does not modify 
the impact of the T-dependence of the drag coeffi¬ 
cient discussed in this letter. 

The results shown have been obtained evaluating 
the drag F from the respective models and then the 
diffusion coefficient D from FDT. Several other op¬ 
tions are possible like evaluating the diffusion from 
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the scattering matrix and the drag from the FDT or 
employing both drag and diffusion from the scatter¬ 
ing matrix. We have seen that while these different 
options may lead to some differences, once they are 
tuned to Raa, the differences in the elliptic flows 
stays within a 10% and in particular our main result 
on the impact of the T dependence of the drag is not 
affected by it. 

In summary, we have evaluated the drag and dif¬ 
fusion coefficients of the heavy quarks within four 
different models. With these transport coefficients 
and heavy quark initial distributions we have solved 
the Langevin equation. The solution of Langevin 
equation has been used to evaluate the nuclear sup¬ 
pression factor, Raa, and elliptic flow, V2 ■ The 
results have been compared with the experimental 
data both at RHIC and LHC energies. Our primary 
intent is to highlight how the temperature depen¬ 
dence of the interaction (drag coefficient) provides 
an essential ingredient for the simultaneous repro¬ 
duction of the nuclear suppression factor, Raa, and 


elliptic flow, V 2 which is a current challenge almost 
for all the existing model. Our work shows that the 
reproduction of the data on Raa(pt) only cannot be 
used to determine the drag coefficient T(T) of heavy 
quarks. We find that the different T-dependences 
of the drag coefficients in the literature can lead to 
differences in V2 by 2-3 times even if the Raa is very 
similar. Our study suggests the correct temperature 
dependence of the drag coefficient cannot be a large 
power of T, like T 2 as in pQCD or AdS/CFT. We 
remind that T(T) nearly constant or weakly decreas¬ 
ing with T would be more typical of a liquid and not 
of a gas. 
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